Rotational excitations in two-color photoassociation 
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We show that it is possible to excite higher rotational states J > 2 in ultracold photoassociation by 
two laser fields. Usually higher J states are suppressed in photoassociation at ultracold temperatures 
in the regime of Wigner threshold laws. We propose a scheme in which one strong laser field drives 
photoassociation transition close to either J — 1 or J = 2 rotational state of a particular vibrational 
level of an electronically excited molecule. The other laser field is tuned near photoassociation 
resonance with J > 2 rotational levels of the same vibrational state. The strong laser field induces 
a strong continuum-bound dipole coupling. The resulting dipole force between two colliding atoms 
modifies the continuum states forming continuum-bound dressed states with a significant component 
of higher partial waves in the continuum configuration. When the second laser is scanned near the 
resonance of the higher J states, these states become populated due to photoassociative transitions 
from the modified continuum. 

PACS numbers: 34.50.Cx, 34.50.Rk, 42.65. Dr, 33.20.Sn 
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1. INTRODUCTION 

Photoassociation (PA) spectroscopy [l|, [|J of ultracold 
atoms by which two colliding atoms absorb a photon to 
form an excited molecular state is an important tool 
for studying ultracold collisional properties at the in- 
terface of atomic and molecular states. PA is particu- 
larly usefuHbr^ producing translationally cold molecules 
;enerating optical Feshbach res- 
1411 - More than a decade ago, 




onance 



and 
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theoretical models [iji [l6j were developed to explain PA 
line shape in the weak-couplin g regim e. The effects of 
laser intensity on PA spectra [13, IS d H HI, [H| 
have been an important current issue. Over the years, 
two-color Raman type PA has emerged as an impor- 
tant method for creating translationally cold molecules 
in the ground electronic configuration. Recently, using 
this method, cold polar molecules [7[ in rovibrational 
ground state have been produced. Molecules created 
by one- or two- color PA of ultracold atoms generally 
possess low-lying rotational levels J < 3. Motivated by 
recent experimental observation of excitation of higher 
rotational states in ultracold PA with an intense laser 
field [23|, we here expore theoretically the possibility 
of rotational excitations in two-color PA. This may be 
important for producing translationally cold molecules 
in selective higher rotational states. Previously, two- 
color PA has been investigated in different other contexts 
EE EE E3, EE ElJpaJuL such as photo-ionization 
of excited molecules [24 l25l l26l. |27| , shielding of atomic 
collision [28l [2Q| . [30| . measurement of s-wave scattering 
length [31(, etc. 

In this paper we propose a method of two-color 
photoassociation of two homonuclear atoms for excit- 
ing higher rotational levels. Our proposed method is 
schematically shown in Fig.l. Laser La is a strong field 
and the laser Lb is a weak one. La is tuned near either 
J A = 1 or J a = 2 rotational state of a particular vibra- 
tional level v of the excited state. This rotational state 



is predominantly accessed by PA transition from s-wave 
scattering state. A photon from La causes PA excita- 
tion from the continnum (s-wave) to the bound level J a- 
A second photon from the same laser can cause a stimu- 
lated de-excitation back to the continnum state. This is a 
stimulated Raman-type process which can lead to signif- 
icant excitation of higher partial waves in the two-atom 
continnum. Now, if a weak laser Lb is tuned near Jb > 2 
states, these higher rotational states get excited due to 
PA from the modified continnum. In this scheme of two- 
color PA, three photons are involved. This does not fit 
into a standard A or V-type process. Here bound-bound 
transition is absent. All the transitions are of continnum- 
bound type. This scheme may be viewed as a combina- 
tion of A and V-type process with continnum acting as an 
intermediate state for V-type transition. In the previous 
Raman-type PA experiments, excited molecular state is 
used as an intermediate state. Furthermore, usually two- 
color PA is carried out in the weak coupling regime. In 
contrast, our proposed scheme involves necessarily one 
strong laser field for inducing strong PA coupling. We 
demonstrate excitation of higher rotational levels in two- 
color ultracold PA by resorting to a simplified model. We 
first evaluate higher partial wave scattering states modi- 
fied due to strong photoassociative coupling [l4| induced 
by the strong laser La- We employ these modified wave 
functions to calculate two-color stimulated line widths 
which are significantly enhanced compared to those in 
the case of one-color. 



The paper is organized as follows. In the following 
section we describe the formulation of the problem and 
its solution. The numerical results and discussion has 
been given in Sec. 3. Finally the paper is concluded in 
Sec.4. 
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r (internuclear distance) 



FIG. 1: A schematic diagram showing the strong (double- 
arrow thick line) and weak (single-arrow thin line) field cou- 
plings between the excited rotational levels J a (a — A, B) and 
the continnum state. Strong laser La modifies the continnum 
state by a two-photon process (curly lines) as described in 
the text. The laser Ls is tuned near resonance with the ro- 
tational levels Jb > 3 which are then populated due to PA 
transition from the modified continnum. Molecular rotational 
levels J = 1 and J — 2 are accessible from s-wave (I = 0) 
scattering state, but J > 3 can only be accessed from higher 
partial-wave (I > 0) scattering states. 



2. THE MODEL AND ITS SOLUTION 

To start with, let us consider that PA laser couples 
continuum (scattering) states of collision energy E = 
ft 2 fc 2 /(2/x) (where /i is the reduced mass) of two alkali- 
type homo-nuclear ground state S atoms to an excited 
diatomic (molecular) bound state which asymptotically 
corresponds to one ground S and another excited P 
atom. Under electric dipole approximation, the inter- 
action Hamiltonian can be expressed as 

H int = E L*-*i 

where di = — er^ is the dipole moment of i-th atom whose 
valence electron's position is given by Vi with respect to 
the center of mass of this atom. Here e represents an 
electron's charge, El is the laser field amplitude and it 
is the polarization vector of the laser. The total Hamil- 
tonian in the center-of-mass frame of the two atoms can 



be written as 

h 2 h 2 

H = H elec (r 1 ,r 2 ;r a , r b ) - — V 2 - ^j^r + H hf + H mt 

(2) 

where H e i ec is the electronic part of the Hamiltonian 
which includes terms corresponding to kinetic energy of 
the two valence electrons, mutual Coulomb interactions 
between nuclei and the electrons, exchange and electronic 
spin-orbit interaction. Here r a and represent the posi- 
tion vectors of the nuclei of atoms a and b, respectively; 
V r and Vr denote the Laplacian operators correspond- 
ing to the relative coordinate r = r a — n and the center- 
of-mass coordinate R = (r + rt)/2 and Hhf stands for 
the hyperfine interaction of two atoms. Under Born- 
Oppenheimer approximation, while solving the electronic 
part of the Hamiltonian, the nuclear coordinates appear 
merely as parameters. PA laser couples only two elec- 
tronic molecular states which are the initial ground and 
the final excited diatomic states represented by (ri, r 2 ; r | 
g) = (/} g (r 1 ,r 2 ;r) and (ri,r 2 ;r | e) = (j} e (r 1 ,r 2 ;r), respec- 
tively. These internal electronic states have parametrical 
dependence on the internuclear coordinate r. They sat- 
isfy the eigenvalue equations 

HeZ ec <Mri,r 2 ;r) = V a (r)<f> a (ri,T 2 ;r); a = g, e 

(3) 

We assume that the matrix element (e | Hi nt \ g) ~ 
fl eg (r) depends only on separation r. Then the center-of- 
mass motion gets decoupled from relative motion. Hence- 
forth we consider only the relative motion. By specifying 
the electronic parts of both the bound and the continuum 
states, one can calculate the matrix element of Hi nt over 
the electronic parts of the two molecular levels involved 
in free-bound transition and thus obtain molecular cou- 
pling strength A(r). 

The continnum-bound dressed state can be written as 
* B (ri,r 2 ;r) = J2a= g ,e®°<( r ) I a ) which is assumed to be 
energy normalized with E being the energy eigen value. 
In the absence of atom-field interaction Hi nt , the prob- 
lem is to find out multi-channel scattering wave-function 
in the ground electronic configuration. The scattering 
channels correspond to two separated atoms a and b in 
hyperfine spin f a and ft,, respectively. The molecular hy- 
perfine state is characterized by the spin F = f a + f b . 
A channel is defined by the angular state | F;f a ,fb,£) 
where ^ = ^ + 7 > = 7 Q + 76 + ~? where £ is the 
mechanical angular momentum of the relative motion of 
the two atoms. This asymptotic basis | T\ f a , fb,fy can 
be expressed in terms of the adiabatic molecular basis 
| T\ S, I, £) [32| , where S and I are the total electronic 
and nuclear spin angular momentum of two atoms. In 
the case of excited molecular state, 5 should be replaced 
by electronic angular momentum J e = S + L. Alterna- 
tively, the adiabatic basis | T; J e ,I,£) can also be ex- 
pressed in terms of | T; J(J e ,£)I). Thus the rotational 
state of a diatom can be expressed in terms of the ma- 
trix element | JttM) = i J J '^rT>^ a (f) where M and 
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n are the z-component of J in the space-fixed and body- 
fixed coordinate frame and f represents the Euler angles 
for transformation from body-fixed to space-fixed frame. 

Djjjj(r) is the rotational matrix element. For ground 
electronic configuration, we have J = £, M = rri£ and 
fl = 0; thereby, fi^M^) reduces to spherical harmon- 
ics Yi mi , We thus express the ground state $g(r) in the 
following form 



t,m t 



(4) 



where ipE'emiir) is the energy-normalized scattering 
state with collision energy E' and j3' E is the density of 
states of unperturbed continnum. Similarly, for a partic- 
ular value of f2, we can expand the excited state $ e ( r ) hi 
the following form 



$ e (r) oc r- 1 [Mr) I JOM}] 



(5) 



M 



Substitution of Eqs. ([!]) and (J5J) into time-independent 
Schrodinger equation leads to coupled differential equa- 
tions. These equations are solved by the use of real space 
Green's function. The detailed method of solution for a 
model problem is given in Appendix A. In our model 
calculations, we consider only a single ground hypcrfinc 
channel. The solution <f) v j(r) can be expressed as 



4>vj(r) = p E , V AjM;t m A vAr)dE' 
Je> w 



(6) 



£m t M 



where <fivj( r ) 1S the excited molecular state (unit- 
normalized) in the absence of laser field and 



A 



J,M;£,m e 



fj,M;£m e 



E s j h ; ft Aj 



H5 + E - E v j + ifry/2 



(7) 



is the probability amplitude of excitation of J from a 
particular partial wave £. Here 



fjM; 



<iP vJ {r')Kj M , tmt (r')^%r')dr' (8) 



is the continnum-bound dipole matrix element and 
AjM-.imt — { JMQ I Q eg | imgO). iPe^ 9 i r ') represents 
the £-th partial wave regular scattering solution in the 
absence of laser field and 

Ef; ft = ttJ J dr'dr^ vJ (r')A JM;ime (r') 

x [Ki(r',r)]A <m(iJM (r)^(r) (9) 

is the partial light shift of the excited state. Here /Q(r, r') 
is the propagator as defined in the Appendix A. The total 



TABLE I: Numerically calculated rotational energies E v j (in 
unit of GHz) and total shift E B ^ (in unit of MHz) for one- 
color laser intensity 7=1 kW/cm 2 for vibrational state v = 
48 of l g excited state. Also given are the rotational energy 
spacings Aj — E v j — E v j-\ (in unit of GHz) for a few lowest 
J values. 



J — ► 


1 


2 


3 


4 


5 


6 


E vJ (GHz) 
Aj (GHz) 
-Ef m (MHz) 


0.57 
19.69 


2.13 
1.56 
22.79 


4.77 
2.64 
17.36 


8.55 
3.78 
11.61 


13.03 
4.48 
6.83 


18.34 
5.31 
2.22 



probability amplitude of excitation Aj for a particular J 
is given by 



.4 



;£.rri£ 

^ h8 + E- e vJ + ihq/2 - Ef lft 



(10) 



i,mi,M 



where Ej = 



E 



E S j£ is the total energy shift 



e,m e ,M J ~'J£ 



of the excited level, fry/ 2 is the natural line width of the 
excited molecular state, E v j is the bound state energy 
corresponding to the bound state solution </>"j of the ex- 
cited state. 8 — lul — ^>A is the frequency off-set between 
the laser frequency lul and atomic resonance frequency 
u>A ■ The ground state scattering solution in the presence 
of PA laser is given by 



1pE£m e (r) = IpsT 9 + ^2 AjM-£'rn e , (E) 

x / K,i{r,r')A tmi ,j M {r')cj>lj{r')dr' (11) 



In the asymptotic limit (r — > oo), the modified scattering 
wave-function behaves like 



Ei'rn, 



cos r}e^m e9 + sin ^^eT ( 12 ) 



where iP'e" 1 ' is the irregular wave function of £-th partial 
wave. Here r\\ is the phase shift due to the applied laser 
field and is given by 



-7T AjM;t'm e i(E)ftme;JM 
jJM;i'm,i 



(13) 



,M 



■pshift 



e,m ll ,M M + E - E ^ + i ^l 2 - E J 
X famt-JM (14) 

where f imt ;jM = ' 4>lj{r')Ai mi .j M (r')^ E ^ eg '(r')dr' . 

(2) 

The two-color partial stimulated line width Tj^ e for a 
particular rotational state Jb is given by 



.(2) 
J B l 



2tt 



(j>° v j{r)AjM;£m e (r)ipE£ me (r)dr 



(15) 



and the total stimulated line width is T 



E 



F (2) 
lm t M 1 J B £ 



(2) 
Jb 



The excitation of particular rotational 
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FIG. 2: Two-color partial stimulated line widths Tf\ (in 
unit of MHz) as a function of AE v .j A (in unit of GHz) at 
collisional energy E = 10 fiK. The intensity I a of laser La 
tuned near Ja = 1 (a) and J a — 2 (b) is 40 kW/cm 2 and the 
intensity Ib of weak laser Lb is 1 W/cm 2 . The total shift 
of the rotational state J a = 1 and J a = 2 are -0.79 



shift 

j Ja 



GHz and -0.91 GHz, respectively. 



state J from the partial wave £ is governed by the fol- 
lowing selection rule 



<£< 



J+\t + ~g\ 



(16) 



where L is the total electronic orbital angular momentum 
L = Z i + ^ 2 and S is the sum of two individual atomic 
spin, i.e. if = "s*i + ~7? 2 - So the lowest possible partial 
wave £ which can make the largest contribution to the 
excitation of rotational state J = 1, 2, 3, 4, 5, 6 are 0, 0, 
1, 2, 3, 4, respectively. The two-color photoassociation 

(2) 

rate K y p ' A for Jb > 2 is defined as 



K 



(2) 



P. I 



(2) 



{VrelCTj) 



hQ T 



npfh- 



0E dE (17) 



7 r ( 7 2) /[(M B +E- E vjB Ef^f + ( 7 + 
?ifc //i is the relative velocity of two 

Jb 



where P Jg 
if } ) 2 /4] and v rel -- 

atoms, a, j — hPj 2 J /k 2 is the inelastic cross-section due 
to loss of atoms. Here (• • •) implies an averaging over the 
distribution of initial velocities, Qt = {^I^KsT/h 2 ) 3 ^ 2 
is the translational partition function and ft = {KgT)^ 1 . 
In the next section, we apply this formalism to a model 
system and obtain numerical results. 



3. RESULTS AND DISCUSSION 

For numerical illustration, we consider a model system 
of two cold ground state (£1/2) 23 Na atoms undergoing 
PA transition from ground state 3 S+ to the vibrational 



TABLE II: Tabulated are one- and two-color partial stimu- 
lated line widths rfg and Fj 2 ^ f at E = 10 /iK for two Sa 
values. Here laser La is tuned near J a — 1 rotational state. 
The intensities of the two lasers are Jb = 1 W/cm 2 and I a = 
40 kW/cm 2 . 





S A = -1.25 GHz 


5 A = -1.48 GHz 


J 


£ 


Tf e (MHz) 


rflt (MHz) 


Ifii (MHz) 


3 


1 


1.55xl0" 05 


0.0158 


0.0107 


3 


2 


0.0000 


0.0111 


0.0065 


3 


3 


0.0000 


0.0085 


0.0051 


4 


2 


0.0000 


0.0128 


0.0075 


4 


3 


0.0000 


0.0072 


0.0047 


5 


3 


0.0000 


0.0103 


0.0061 




20 25 

r (a.u.) 



20 25 30 

r (a.u.) 



FIG. 3: The subplots (a), (b) and (c) show the light-induced 
scattering wavefunctions ipEim e (in unit of Bohr radius -1 ^ 2 
x Hartree" 1/2 ) for p- {£ = !), d- {£ = 2) and f-wave (1 = 3), 
respectively. The solid and dashed curves correspond to the 
detunning 8a ~ —1.25 GHz and 5a = —1.48 GHz, respec- 
tively. The subplots (d), (e) and (f) exhibit the corresponding 
field-free regular wavefucntions i^^ 9 ■ All the wave functions 
are plotted at collisional energy of 10 /iK and intensity Ia = 
40 kW/cm 2 . 



state v = 48 of the excited molecular l g state [22 



At 



large internuclear distance this l g potential correlates to 
2 Si/ 2 + 2 P%/2 fr ee atoms and at short r ang e to 1 1 n g 
Born-Oppenheimer potential. In Ref. [23j higher ro- 
tational lines upto J = 6 have been clearly observed 
in PA with an intense laser field. The centrifugal bar- 
rier of £ > of the two- atoms lies at r > 50ao (ao= 
Bohr radius) whereas PA excitations occur at r ~ 27eto. 
Therefore, the higher rotational states will be unlikely to 
be populated by PA transitions from £ > partial-wave 
scattering states at ultra-cold temperatures in the weak- 
coupling regime. Previously, higher rotational levels have 
been excited in PA spectroscopy due to resonant dipole- 
dipole interaction with transition occurring at large in- 
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FIG. 4: tanr/f (rjf is the light induced phase-shift) is plotted 
as a function of detunning AE V] j A =2 (MHz) when La is tuned 
near J A = 2. The total shift Ef^Ll at 40 kW/cm 2 is -0.91 
GHz. The other parameters are I a = 40 kW/cm 2 and E — 
10 fjK 



ternuclear separations [33|, [3J]. The numerically calcu- 
lated rotational energies E v j, energy shifts Ej** 1 and 
the corresponding energy difference Aj = E v j — E v j_\ 
for six lowest J values are given in Table I. To demon- 
strate the working of our proposed scheme, we resort 
to a simplified two-state calculation. We consider only 
one ground hyperfine channel with F — 4, f a — 2 and 
fb = 2 in the absence of any external magnetic field. In 
the excited molecular state, we neglect the hyperfine in- 
teraction. The two-color partial stimulated line width 

(2) 

is plotted as a function of detunning AE v> j A = 

hS A + E- E Vt j A - Ej h A ift in Fig.2 for J B ranging from 3 
to 6. The strong laser La is tuned near J a = 1 (Fig.2a) 
and J A = 2 (Fig.2b). From Fig.2 we notice that rf^ t 
strongly depends on the detunning AE v- j A of the strong 
laser from PA resonance of the rotational level J a- The 

(2) 

maximun of Lj^ e occurs at AE v j A = 0. For lower Jb 
values, the probability of rotational excitation is higher. 

For comparison, we also calculate one-color partial 



stimulated line widths T 



(0) 



sion r 



jg for J > 2 from the expres- 
^\fjM;im e \ 2 ■ The one-color total stimulated 
line width is Ty = m m^j}- ^ ^ energy and 
at laser intensity 1 W/cm 2 , the one-color partial stim- 



(0) 
.]( 



ulated line widths T 



(0) 
J=3f=l 



15.46 Hz, T 



(0) 

J=Al=2 



0, 



r jL'a f=3 — 0. A comparison between one- and two-color 
partial stimulated line widths has been made in Table II 
for Jb > 2 at collisional energy 10 /iK. The two-color 
total line widths T^ =3 , r^ =4 , r^ =5 are 0.03537 MHz, 
0.0200 MHz and 0.0103 MHz, respectively when 5 A = 
-1.25 GHz and they are 0.02229 MHz, 0.0122 MHz and 
0.0061 MHz, respectively for 5 A = -1.48 GHz. The corre- 
sponding one-color weak-coupling partial as well as total 




4 6 
E (MHz) 



FIG. 5: The Two-color total stimulated line width (in 

J B \ 

unit of MHz) for different Jb (as indicated in the plots) is 
plotted as a function of collisional energy E (in unit of MHz) 
when La is tuned near J A — 1 for 8a = -1.25 GHz (a), 8 A = 
-1.48 GHz (b) with I A = 40 kW/cm 2 and I B = 1 W/cm 2 . 



TABLE III: Tabulated are the tan% when the laser La is 
tuned near J A — 2 at E — 10 /iK for three values of 5a- The 
parameters are Ia = 40 kW/cm 2 , E v j A =2 = —2.138 GHz 
0.91 GHz. In the field-free case, tan?7° =1 

i and 



and Effll 



= -1.57x10" 
t&nr]i =4 ~ 



tan rfl- 



= 1.20x10" 



tan rfl- 





8 A = -2.95 GHz 


8 A = -3.049 GHz 


8 A = -3.17 GHz 


i 


taxirjf x 10 4 


tan??]' x 10 4 


tanjjf x 10 4 


1 


263.00 


36900.00 


-229.00 


2 


4.93 


659.00 


-4.09 


3 


0.01 


2.01 


-0.01 


4 


0.00 


0.01 


0.00 



stimulated line widths T ?} and TK'' for the same rota- 



(0) 



Jt 



J 



tional states with laser intensity of 1 W/cm are van- 

(2) 

ishingly small while the two-color partial L and total 
Tj 2 ^ exceed and by several orders of magnitude. 

We find the energy shift \Ej is 0.79 GHz which ex- 
ceeds the spontaneous line width 7 (say 2 MHz for model 
calculation) by two orders of magnitudes. 

In order to trace the origin of increment of Ij^ we plot 
perturbed ipEem e for i ^ when laser La is tuned near 
J A = 1 and the corresponding field-free regular func- 
tions V-'bJ 69 i n Fig.3. It is clear from this figure that 
the amplitudes of i\) Eime are greatly enhanced by several 



orders of magnitude than that of tp 



O.reg 
Ei 



Next, we cal- 



culate tan rjf by using Eq. (jT4j) when laser La is tuned 
near J a — 2. These are given in Table III for 6 a = - 
2.95 GHz, -3.17 GHz and -3.049 GHz. The first two 5 A 
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-40 -30 -20 -10 10 20 
AE vJ (MHz) 



4. CONCLUSION 

In the present paper we have developed a two-color 
PA scheme for the excitations of higher (Jb > 2) ro- 
tational levels which are generally suppressed in the 
Wigner threshold law regime. We have calculated two- 
color stimulated line width (for Jb > 2) by fixing strong 
laser either near J a — 1 or J a = 2 state and tuning 
another weak laser to higher rotational ( Jb — 3, 4, 5) 
states. Then we have compared these with one-color line 
widths. The enhancement of stimulated line width is a 
result of strong-coupling photoassociative dipole interac- 
tion which in turn modifies the continuum states. This 
proposed method may be important for coherent control 
of rotational excitations and manipulation of optical Fes- 
hbach resonance of higher partial waves. 



APPENDIX A 



FIG. 6: The upper panel (a) shows two-color photoassociation 
rate Kp A (in unit of meter 3 sec -1 ) as a function of atom-field 
detunning 5b (in unit of GHz) for three higher rotational 
levels Jb (as indicated in the plots) for 5a = -1.25 GHz when 
the laser La is tuned near J a = L The lower panel (b) 
shows the same but as a function of detunning AE v ,j B — 
hS B -E v ,j B -Ej h J ft (in unit of MHz) from PA resonance. The 
other parameters for both the panels are I a = 40 kW/cm 2 , 
Ib = 1 W/cm 2 and T = 100 /iK. 



The mathematical treatment given here is closely re- 
lated to our earlier work [TJ]. Treating the laser field 
classically, the effective interaction Hamiltonian under 
rotating wave approximation the two state basis can be 
expressed as 



H l e f f = e X p(-iSt)Q eg (r)\e)(g\ +H, 



(Al) 



Form time- independent Schrodinger equation H^e = 
E^Sf e, we obtain two coupled equations 



values correspond to off-resonant and the last one to res- 
onant condition. The variation of tan 77^ with AE v j A= 2 
is plotted in Fig. 4 which exhibits resonance for higher 
partial waves induced by strong-coupling PA. The en- 
hancement of the partial (I ^ 0) wave amplitude is due 
to the term Y,t>m t ,M A JM;l' mi , (E) of Eq. (TTT]). In fig- 

(2) 

ure 5, the two-color total stimulated line width Tjj is 
plotted as a function of collisional energy E for two off- 
resonant 8 a values when L^ is tuned near J a = 1- The 
magnitude of r^j for higher rotational states (Jb — 4, 5) 
is less than that of Jb = 3. This is due to the fact that 
the lowest possible partial wave contribution to the ex- 
citation of rotational states Jb = 4 and Jb — 5 are d 
and /, respectively while Jb = 3 state can be populated 
from p wave which has rotational barrier lower than that 
of d- and f-wave. The two-color photoassociation rate 
K pA as defined in Eq. (frT|) has been plotted as a func- 
tion of 5b (Fig. 6a) and AE Vi j b (Fig. 6b). The spectra 
in Fig. 6b are red-shifted due to the presence of the term 
Tj^ in Eq. (JT7J). From the selection rule, it is obvious 
that Jb — 3,4,5 rotational states can not be populated 
by a PA transition from s-wave scattering state. But the 
appearance of the Jb — 3,4,5 lines in PA spectra is an 
indication of the significant modification of the partial 
scattering wavefunctions by intense light field. 



l-^ 2 r + V g (r)-E]%(r) 



{-^V 2 r + V cx (r) - E - h6]t> e (r) 



-0, e (r)$ e (r) (A2) 



-ft eg (r)^ g (v)(A3) 



Here V g is assumed to include hyperfine interaction of the 
chosen channel. Substituting Eqs. ([!]) and (f5]) into the 
Schodinger equations ()A2j) and (|A3|) we get two coupled 
equations 



+ Bjij) + Vcx{r) ~h5-E~ih± 



4>vj 



2fi dr 2 



+ B e (r)+V s (r)-E 



M 



A(.m t -JM<t)vJ 



(A5) 



where Bj(r) = ti 2 / (2[ir 2 )[J( J+ 1) - ft 2 ] is the rotational 
term of excited molecular bound state in the absence 
of nuclear spin, Bt(r) — ti 2 / (2/ J ir' 2 )£(£ + 1) is the cen- 
trifugal term in collision of two ground state (S) atoms, 



7 



■^Etm t { r ) = f E > 0E>ipE>tmt(,r)dE' . The above two equa- 
tions are solved by the green's function method by setting 
0. The single channel scattering 



A 



= A 



tmi;JM 



equation becomes 



■^ + B i (r) + V t (r)-E 



= 0. (A6) 



Let 4>Ei e9 (r) and ^gj~ r (r) represent the regular and ir- 
regular solutions of the above equation. The appropriate 
Green's function for the scattering wave function can be 
written as 

Ki(r,r') = -TV^% e9 {r)xl)% l [ r {r') {r > r) (A7) 



Using this Green's function, we can write down the solu- 
tion of equation (|A4p in the form 



where 



4>vj[r) = - dr'A JM ;tmi{r')G v (r,r')iij£ me (r') 



f /3 E ,J2 A JM;im e <t>vAr)dE' (A13) 



K t {r, r') = -^% e9 {r')^[ r {r) (r' < r) (A8) 

The regular function, ^eT 3 ( r ) vanishes at r=0 and the 
irregular solution iPeT P ( r ) * s defined by boundary only 
at r — > oo. The energy normalised asymptotic form of 
both regular and irregular wave function is 



AjM-,em e = ^2 j dr' KjM;im l {r')4>ij{r')'ilJEim e {r') 

X h5 + E-E vJ + ih 1 /2 ( ' A14 ^ 
Substituting equation (|A13|) into equation (|A5|) we ob- 
tain 



/0,reg _ 

Vei — 



2/i . ., in n 
__ sln(fcr __ + ^) ) 



oo (A9) 



0,irr 
El 



2n, £tt 
^^cos(fcr- — +4), r^oo(AlO) 



IpEtme (r) 



E A JM;e'm e , A trnr.JM(r)(f>vAr) 
I'm. i M 



(A15) 



where rfl is the phase-shift of £-th partial wave in the 
absence of PA coupling. The homogeneous part of (|A4|) 
with 7 = is 



The scattering solution can now be expressed as 



b° vJ = {h8+E)4>lj = E vJ <f>°j *l>Et mt (r) = 4>e7 9 + A JM;e'm e ,(E) 



(AH) 

The Green function corresponding to these rovibrational 



states (f>® j can be written as 



x / JC e (r,r')A eme .j M (r')<P v j(r')dr(Al6) 



G v (r,r') 



1 



H6 + E — E vJ + ihj/2 



(A12) 



On substitution of equation (|A16[) into (|A14[) and after 
some algebra, we obtain 



iJ,M;l,m e 



hS + E- E vJ + itvy/2 - E S J 



shift 



JM-.tmt 



E 



shift 
J£ 



E 

t'^l,m i 



A 



JM'-.l "m. 



(A17) 



M' 



Let D = HS + E - E vJ - E s ]f St + ifvy/2. Now, adding 
a term D^ 1 E S J^ Aj y M-,e,mi on both side of equation 
(|A17[) . we can express Aj^M-e.me in terms a quantity 
Aj = J2e me M Aj,M-e,m e as well as other parameters. On 



summing over all possible £, me, M we can evaluate Aj. 
Having done all these algebra, we can explicitly express 

fj,M;tm t + Ej^ 1 1 Aj 
Aj t M;l,mt = M + E _ EvJ+lhl/2 ( A18 ) 
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and 



fjM;lm e 



shift 



(A19) 
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